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ABSTRACT 

We investigate the dynamics of an Alfven surface (where the Alfven speed equals the advection 
velocity) in the context of core collapse supernovae during the phase of accretion on the proto-neutron 
star. Such a surface should exist even for weak magnetic fields because the advection velocity decreases 
to zero at the center of the collapsing core. In this decelerated flow, Alfven waves created by the 
standing accretion shock instability (SASI) or convection accumulate and amplify while approaching 
the Alfven surface. 

We study this amplification using one dimensional MHD simulations with explicit physical dissi- 
pation (resistivity and viscosity). In the linear regime, the amplification continues until the Alfven 
wavelength becomes as small as the dissipative scale. A pressure feedback that increases the pres- 
sure in the upstream flow is created via a non linear coupling. We derive analytic formulae for the 
maximum amplification and the non linear coupling and check them with numerical simulations to 
a very good accuracy. Interestingly, these quantities diverge if the dissipation is decreased to zero, 
scaling like the square root of the Reynolds number, suggesting large effects in weakly dissipative 
flows. We also characterize the non linear saturation of this amplification when compression effects 
become important, leading to either a change of the velocity gradient, or a steepening of the Alfven 
wave. 

Applying these results to core collapse supernovae shows that the amplification can be fast enough 
to affect the dynamics, if the magnetic field is strong enough for the Alfven surface to lie in the 
region of strong velocity gradient just above the neutrinosphere. This requires the presence of a 
strong magnetic field in the progenitor star, which would correspond to the formation of a magnetar 
under the assumption of magnetic flux conservation. An extrapolation of our analytic formula (taking 
into account the nonlinear saturation) suggests that the Alfven wave could reach an amplitude of 
B ~ 10^^ G, and that the pressure feedback could significantly contribute to the pressure below the 
shock. 

Subject headings: waves — accretion — magnetic fields — magnetohydrodynamics — supernovae: general 



1. INTRODUCTION 

The theoretical mechanism responsible for the explo- 
sion of massive stars has become a multidimensional 
puzzle when spherical model s were ruled out by ac - 
curate numerical simulations (jLiebendorfer et aD 120011 ) . 
In all mechanisms proposed recently, multidimensional 
fiuid motions are essential for s uccessful explosions of 
massive stars via core co llapse (|Marek fc Jankal 120091 : 
iBurrows et al.l 120061 12007D . These motions are induced 
by several hydrodyn amic instabilities: c onvection in the 
proto-neutron star (jPessart et al.l [20061 ) . convection in 
the neutrino hea ted region between the neutrinosphere 
and the shock (jHerant et al.l I1992D and the standing 
accret i on shock instability (SASI, e.g. iBlondin et al.l 
(I200I : iFodizzo et al.l (l2007h l 

The effect of a magnetic field on the dynamics has 
been studied mainly in the context of very rapid ro- 
tation. In this case the magnetic field can be ampli- 
fied to a dynamically significant strength by the winding 
of the initial field and by the ma gneto-r otational insta - 
bility (jAkivama et al.l 120031 : lObergaulinger et al.l I2OO90 . 
This magnetic field could then extract enough rota- 
tiona l energy to power an explosion driven by magnetic 
jets (jBisnovatvi-Kogan et aLlUQTet iShibata et all 120061 : 

|jerome.guilet@cea.fr| 



IBurrows et"alll2007[ ). 

For moderate rotation rates, the magnetic field is usu- 
ally neglected because the magnetic pressure is thought 
to be much smaller than the thermal pressure. The cri- 
terion determining the significance of the magnetic field 
effects, however, m ay be more complex. For example, 
IGuilet fc Foglizzd (|2010) suggested that a magnetic field 
can affect SASI significantly if the Alfven speed is com- 
parable with the advection velocity, which is a much 
less stringent criterion in the very subsonic part of the 
flow. This may explain the surprising stabilization of 
SASI by a radial magnetic field which pressure is negli- 
gi bly small {P^^./P < IQ-^), in t he model 2DB14Ain 
of lEndeve et all poTol ). Moreover. lEndeve et al.l (poTof ) 
showed that even in the absence of rotation, the magnetic 
field could be amplified by the turbulent flow driven by 
SASI. This motivates our study of the effect of a mod- 
erate magnetic field in the context of negligible rotation, 
when the Alfven and the fluid velocities are comparable. 

The Alfven surface is defined as the surface where the 
flow velocity equals the Alfven velocity. Such a surface 
exists in the context of a wind or a jet ejected from a mag- 
netized object (e.g. the Sun) because the fluid accelerates 
from a subAlfvenic to a superAlfvenic speed. Symmet- 
rically in the context of accretion of a magnetized fluid 
onto a solid object, if the flow is superAlfvenic far from 
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the accretor it has to decelerate to a subAlfvenic speed 
before setthng down onto the surface. The fundamental 
difference between these two situations is that an Alfven 
wave that is propagating against the flow accumulates at 
the Alfven surface in a decelerated flow (as its velocity is 
decreasing to zero), while it diverges from it in an accel- 
erated flow. This is similar to the hydrodynamical situa- 
tion at the sonic point: in a decelerated transsonic flow, 
the accumulation of acoustic waves leads to the formation 
of a shock wave, while no such thing happens in an ac- 
celerated flow. The formation of a transAlfv enic discon- 
tinuity was shown to be non evolutionary by iSvrovatskiil 
(flQSl : in the presence of an arbitrarily small perturba- 
tion, it instantaneously disintegrates into several other 
discontinuities. What happens then at the Alfven sur- 
face? To our knowledge, the only stu dy dealing with a 
decelerated transAlfvenic flow is that of lWilliams (|1975| ). 
He showed that Alfven waves not only accumulate but 
are actually amplified while approaching the Alfven sur- 
face in the sense that their energy flux increases. He 
proposed that this could cause a broad turbulent zone 
below the Alfven surface. 

Such an Alfven surface is expected to exist in the con- 
text of core collapse supernovae during the phase of ac- 
cretion onto the proto-neutron star. This holds even for 
a very weak magnetic field (and thus a very slow Alfven 
speed) because the infall velocity vanishes at the center 
of the sta r. Furthermore strong A lfven waves are created 
by SASl CGuilet k Foglizzoll2oioD and bv th e convection 
inside the proto-neutron star (jSuzuki et al.l 12001) . The 
primary goal of this paper is to study the evolution of 
these waves properties as they reach the Alfven surface 
in the conditio ns prevailing duri ng the collapse of a mas- 
sive star core. iWilliamsl (|f 9751 ) made a linear analysis 
of high frequency Alfven waves in a cold supersonic and 
inviscid flow. We extend his analysis by including the ef- 
fect of thermal pressure, viscosity, and resistivity, and by 
considering also the low frequency regime. Furthermore, 
we focus on the non linear dynamics in a subsonic flow, 
and show the existence of a non linear pressure feedback. 
For this purpose we use a combination of MHD simula- 
tions and analytic arguments, which we apply to a model 
of Alfven surface that is voluntarily as simplified as pos- 
sible. 

In Section 2 we describe our set up and our numeri- 
cal method. Section 3 gives a qualitative description of 
the simulations, while Sections 4 and 5 are devoted to 
a quantitative description of the amplification of Alfven 
waves in two different frequency regimes. The satura- 
tion of the Alfven amplification is studied in Section 6. 
Finally we discuss the astrophysical consequences of our 
results in Section 7 and conclude in Section 8. 

2. THE MODEL 

2.1. Stationary flow 

In order to study the dynamics of an Alfven surface in 
a general perspective, we choose to simplify the physical 
set up as much as possible. As the phenomena discussed 
below take place in the vicinity of the Alfven surface, 
the convergence of the flow is not expected to play a 
crucial role, hence we restrict our study to a planar ge- 
ometry. The perfect gas is described by an adiabatic 
index 7 = 4/3 characterizing a gas of relativistic degen- 
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Fig. 1. — Schematic view of the setup. An external potential 
step located in — Lv < ^ < Lxj decelerates the flow from a super- 
Alfvenic to a subAlfvenic velocity. An Alfven wave propagating 
against the flow is sent from either the lower or the upper bound- 
ary, and accumulates at the Alfven surface. The upper part of the 
figure represents the propagation speed of this wave {v+vj^), which 
vanishes at the Alfven surface located in a:: = 0. On the contrary, 
Alfven waves propagating in the same direction as the flow are not 
affected by the Alfven surface as their propagation velocity t) — da 
does not vanish. 

crate fermions. Cooling processes are neglected for the 
sake of simplicity. The Alfven surface in the stationary 
flow is created through the adiabatic deceleration of a 
superAlfvenic flow to a subAlfvenic velocity by an ex- 
ternal potential step (Figure [T]). The superAlfvenic flow 
located at a; > Ly has a velocity and a magnetic field 
that are aligned along the x direction {v = < 0). The 
external potential is a linear step extending in the region 
— Lv < X < L\/ and with a depth chosen in such a way 
that the Alfven surface is located at a; = 0. The incident 
flow at a; > Lv is denoted with the subscript "in" (e.g. 
the advection velocity Vim sound speed Qn, and Alfven 
speed WAin), the downstream flow at a; < — Lv with the 
subscript "out", and the Alfven surface at a: = with 
the subscript 0. 

The stationary flow profile is determined by the con- 
servation of energy, mass, magnetic flux and entropy: 
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dx 
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where the dimensionless entropy S is defined by: 



S . 



1 



7-1 



log 



P 



(1) 

(2) 
(3) 
(4) 

(5) 



As a consequence, the sound speed, density and Alfven 
speed (defined by va = B j ,JJIqp) scale with the advec- 
tion speed as: 
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The Alfven velocity decreases less than the advection 
velocity, which allows the superAlfvenic fluid to become 
subAlfvenic when compressed by the external potential. 
The potential step has the following depth: 
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(9) 



The free parameters of the stationary flow are thus 
Lv, vo, Win and cq. Other quantities can be obtained 
from these (e.g. VAm = V^ini'o)- 

2.2. Alfven waves 

From the upper or the lower boundary, we send either 
a circularly or a linearly polarized Alfven wave that is 
propagating against the flow. The transverse magnetic 
field and velocity of a circularly polarized wave are de- 
scribed by: 

(5B = SBin [cos {kxX — Lot) Uy + sin {k^x — ujt) u^] ,(10) 



SB 



(11) 



where Uy and are unit vectors along the y and z di- 
rections, 6Bia is the amplitude of the Alfven wave (it is 
noted 5i?out when the Alfven wave is sent from below), 
(jj its frequency and its wave number: 



kx — 
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In order to generate a symmetric displacement with re- 
spect to the initial position, the circularly polarized wave 
is replaced by a linearly polarized one during the first 
quarter of a period {0 < t < T^if/i). A linearly polarized 
Alfven wave is obtained from the same set of equations 
by retaining only the y component of the magnetic field 
in Equation (ITO|) . 

The circularly polarized Alfven wave is an exact so- 
lution of the MHD equations in a uniform flow, in the 
absence of dissipation, if its amplitude SB is constant in 
time and space. However at the time we start to send the 
Alfven wave the transition from a zero to a non zero am- 
plitude is not an exact solution, and would create sound 
waves. To avoid this, we also perturb the density, pres- 
sure, and speed of the stationary flow using a second 
order description of the Alfven wave as is explained in 
Appendix \^ 

2.3. Parameters and dimensionless numbers 

Unless otherwise specified, the physical quantities are 
normalized by setting vao = —vq — ^, Po = ^ and 
= 1- The free parameters of our setup are then: 
Vim Cin, the resistivity 77, and the viscosity v. Alfvenic 
perturbations are characterized by their frequency w and 
their amplitude SBin (or SBout if sent from below). We 
explored the parameter space by varying the parameters 
one by one, using the following baseline values: Cin = 2wo, 
Win = 1.5wo, V = 2.10-^ v = 0, SBin = lO-'^Bo, and 
w = 1. 

We define below a few physical quantities which are 
important for the dynamics of the Alfven surface and 
will be useful in the rest of the paper. The strength of 
the gradients is described by: 



d{v + Va) VAin + ■ 



dx 



(13) 



It is also a typical growth rate of the Alfven wave ampli- 
tude (Section 4.1). Its ratio with the wave frequency de- 
fines a dimensionless number e that determines whether 
the wave is in the high frequency regime (e ^ 1, stud- 
ied in Section 4) or in the low frequency regime (e 3> 1, 
Section 5): 

The illustrative simulation described in Section 3 corre- 
sponds to e ~ 0.15, and has a behavior typical of the 
high frequency regime. 

The ratio of kinematic viscosity ly to resistivity r/ (the 
magnetic Prandtl number) was varied from to infinity, 
and the same results were found in these two limiting 
cases. For simplicity, in most of the paper, we consider 
only the effect of resistivity ; however the same results 
would be obtained with viscosity instead. As shown in 
Sections 4 and 5, a typical length scale for the dissipation 
can be defined as: 



ri + V 



(15) 




Its ratio with the lengthscale of the gradients is related 
to the Reynolds number: 



(16) 



In the simulations presented here the dissipation is varied 
in the range Re ~ 10^ - 10^ (with ~ 10^ for our 
baseline parameters). The core collapse supernovae fluid 
can be rather viscous inside the neutrinosphere due to 
neutrino viscosity and is expected to lie in this range of 
Reynolds number (Section 7). Most other astrophysical 
fluids would be expected to be much less dissipative. 

2.4. Numerical method 



We use the c ode RAMSES (iTevssied [200I 
iFromang et al.l [20061 ) to solve the usual MHD equations 
in one dimension along the x-direction. Note however 
that the velocity and magnetic field are three dimen- 
sional vectors with non zero components along the y and 
z directions. RAMSES is a second order Godunov type 
code that uses the MUSCL-Hancok scheme to evolve the 
MHD equations. For the calculations presented in this 
paper, we used the MinMo d slope limiter alon g with 
the HLLD Riemann solver (jMivoshi fc Kusanol |2005) . 
The computational domain is a one dimensional box 
of length 4Lv, comprised between — 1.5Lv and 2.5Lv 
when the Alfven wave is sent from above, and between 
— 2.5Lv and 1.5Lv when the Alfven wave is sent from 
below. 

After the stationary flow has settled to its numerical 
equilibrium, we start sending an Alfven wave from ei- 
ther the upper or the lower boundary. For this purpose 
we impose in the ghost cells the conditions described by 
Equations (HUll-dni) and (|X3)) - (|M)) . We use a zero gra- 
dient condition for the other boundary, and checked that 
the reflections caused by such a choice are negligible (see 
Section 3). 

The dissipation is taken into account explicitly with a 
kinematic viscosity v and a resistivity rj. For given dis- 
sipation coefficients the resolution is chosen in the range 
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Fig. 2. — Magnetic field lines projected in the plane (xy) at 
the end of the simulation with our baseline parameters (a small 
amplitude circularly polarized incident wave in the WKB regime, 
ui-^ /u) ~ 0.15). The scale in the y direction has been normalized by 
the displacement of the incident wave Syi^ = 5Bin/{kinBo). The 
a;— axis is normalized by the scale of the stationary gradients. 

Ax/Lv ^ 5.10^'' — 4.10^^ SO that the numerical dissi- 
pation is significantly smaller than the physical one and 
does not affect the result: this typically corresponds to 
Ax < Xdis/5. 

3. QUALITATIVE DESCRIPTION OF THE SIMULATIONS 

For the purpose of illustration we first focus on our 
baseline set of parameters and on the case where the 
Alfven wave is circularly polarized and sent from above. 
The frequency such that wy/w ~ 0.15 makes it represen- 
tative of the high frequency regime studied more quan- 
titatively in Section 4. 

As it approaches the Alfven surface, the Alfven wave is 
slowed down to a velocity approaching zero, which causes 
its wavelength to decrease to zero (Figure [5]). At first, its 
transverse magnetic field and velocity are amplified to in- 
creasingly large values while the transverse displacement 
stays roughly constant (Figures [2] and [3|). The Alfven 
wave is efficiently damped when the wavelength becomes 
as small as the dissipative scale, either resistive or vis- 
cous. The resulting amplitude profile of the transverse 
magnetic field peaks close to the Alfven surface, at a 
maximum value which we refer to as (5i?max) and then 
quickly decreases to zero (Figure [3]) . This enables the 
establishment of a stationary state, where the amplitude 
of the Alfven wave is constant in time (but SByz and 
Svyz are oscillating as the wave is propagating). The 
dissipation of the Alfven wave is measured at the lower 
boundary by the change of entropy SS as compared to 
the unperturbed flow (Figure[3l bottom panel). This am- 
pliflcation and dissipation are explained quantitatively in 
Section 4 for the WKB regime, and in Section 5 for the 
low frequency regime. 

The total pressure profile in Figure [3] can be explained 
by two distinct phenomena. The negative extremum (sig- 
nifying a decrease as compared to the stationary pressure 
profile) for a: > is the contribution from the Alfven 
wave, which amplitude has a similar profile. The sec- 
ond order expansion of the Alfven wave explained in Ap- 
pendix gives a good description of the contribution of 
the Alfven wave to the total pressure (dashed line in the 
middle panel of Figure |3]). The second effect is the ere- 
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Fig. 3. — Profiles in the stationary state at the end of the simu- 
lation with our baseline parameters. The stationary profiles have 
been subtracted in order to show the effect of the Alfven wave on 
the flow. Upper panel: y component of the magnetic fleld (full 

line) and transverse magnetic field Btrans = \J + (dashed 

line) normalized by the amplitude of the incident Alfven wave. 
Middle panel: total pressure defined as the sum of the thermal 
pressure and the magnetic pressure Ptot — -P+-B^/(2/io) (full line). 
The contribution of the Alfve n wave to the total pressure has been 
estimated using Eauation lA3l (dashed line). Bottom panel: entropy 
profile. The pressure and entropy have been normalized by SB?^ 
and SBf^/Bg in order to define dimensionless quantities that are 
independent of the incident wave amplitude in the linear regime 
(see Section 6). 

ation of a pressure feedback propagating up and down, 
which can be visualized by the difference between the 
dashed line and the full line. The amplitude of the feed- 
back sent upward (respectively downward) is measured 
at the upper (resp. lower) boundary by (JPgup, Svsup and 
<5psup (resp. SPini, Svini and Sp-mi)- For both signals, 
we verified that the relation between the pressure, ve- 
locity and density perturbations correspond to what is 
expected from an acoustic wave propagating either up or 
down ('Equations [13^5)) . This gives us confidence that 
the boundary condition we apply does not cause any spu- 
rious reflections. We find that the pressure feedback in- 
creases the pressure and decelerates the flow upstream 
of the Alfven surface, while it decreases the pressure and 
decelerates the flow downstream. This pressure feedback 
is explained in Section 4.4. 

The time evolution of the maximum transverse mag- 
netic field, the pressure feedback, and the entropy cre- 
ation (shown in Figure S]) gives some useful information 
on the origin of the pressure feedback. The maximum 
transverse magnetic field increases from the time the 
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Fig. 4. — Time evolution in our baseline simulation of the max- 



imum transverse magnetic field SBy 



(full line), the upstream 



pressure feedback 5Psup measured at the upper boundary (dashed 
line) and the created entropy SS measured at the lower bound- 
ary (dotted line). S-B^niaxi '^^'sup and 5S have been normalized by 
their final value in the stationary state. The time is normalized 
using ojv , a typical growth rate for Alfven wave amplification (see 
Section 4). 

Alfven wave enters the gradients at t ~ 0.5 — 1, un- 
til t ~ 4 when dissipation becomes important. At this 
time the entropy starts to increase and quickly reaches 
its maximum valu43- On the other hand the pressure 
feedback increases soon after the beginning of the ampli- 
fication. This is significantly before the entropy creation 
starts, indicating that the primary source of this pressure 
feedback lies in the amplification of the Alfven wave by 
the gradients rather than in its dissipation. 

4. THE WKB REGIME 

4.1. Non dissipative waves 

When the Alfven wave has a much shorter wavelength 
than the typical scale of the gradients, its evolution can 
be described in the WKB formalism. Close to the Alfven 
surface, this approximation is valid if the parameter e de- 
fined in Eq. (|14|) satisfies e ^ 1/4, which corresponds to 
the high frequency regime. Note that, in the vicinity of 
the Alfven surface, this parameter e does not depend on 
the distance from the Alfven surface although this may 
seem counterint uitive as the wavelength decreases to zero 
(|Williamsl [19751) . This can be understood by the fact 
that the lengthscale of variation of its effective propaga- 
tion velocity also scales like the distance x to the Alfven 
surface. 

In the WKB approximation and in the absence of dis- 
sipation, the wave action of an Alfven wave is conserved 
(heuristically this conservation is equiva l ent to the con- 
servation of wave quanta, see IJacauesI (|1977f )). This 
property can be used to describe the amplification of 
the Alfven wave, thus avoiding the lengthy calculations 
needed for a full WKB approach. The wave action den- 
sity is defined by: 

A^-, (17) 

UJQ 

^ The propagation of entropy and pressure to the boundaries 
where they are measured induces a delay, which is rather short in 
da;/(c + v) O.duj^^ for the pressure 
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dx/v ~ 0.3aj„ for the entropy. 



where E is the energy density and ojq — lu 
frequency in the frame of the fluid: 
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k.v is the 

(18) 
(19) 



The conservation of the wave action can then be stated 
as: 

dA 



dt 



+ V. [(v + wp,)A] = 0. 



(20) 



In a ID stationary state, the wave action flux is thus 
constant: 

{v + va)A = — = constant. (21) 



UJVA 



Mo 



The following scaling of the Alfven wave amplitude is 
obtained: 



SB cx 



VA + V 



(22) 



In the absence of dissipation, the amplitude of an in- 
cident Alfven wave thus increases to infinity when ap- 
proaching the Alfven surface. Such a divergence could 
be expected from a mere accumulation of the wave en- 
ergy: assuming a constant energy flux and a vanishing 
propagation speed implies that the amplitude diverges. 
What happens to the energy flux of the Alfven wave? 
The relation between the energy density and the action 
density (Eq. [T7)) shows that as ujq is varied, both quan- 
tities cannot be conserved at the same time. When an 
Alfven wave approaches an Alfven surface, the increase 
of Wo impl ies that the ene rgy flux of the Alfven wave also 
mcreases (|Wilhamslll9"75[) : 



F = (vA 



SB^ VA 
V) cx . 

Mo VA + v 



(23) 



This shows that the Alfven wave not only accumulates 
but actually amplifies when approaching the Alfven sur- 
face. A wave vector can be associated with this amplifl- 
cation of the energy flux: 



^ _ 1 dF _ I dvA 

amp p 



1 d{v + VA) 



V + VA 



dx 



(24) 



Close to the Alfven surface w-f wa ^ wa, hence the second 
term is dominant and may be approximated by: 



Wv 

V + V\ 



1 

X 



(25) 



The growth rate that can be associated with the Alfven 
wave amplification is: 



^amp 



{vA + w) ~ - 



d{v + VA) 

dx 



(26) 



Thus the steeper the gradient, the faster the ampliflca- 
tion. Interestingly, in the idealized case of a compact 
step, the growth rate would be inflnite. 
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4.2. Dissipation of Alfven waves by viscosity or 
resistivity 

In a uniform flow, an Alfven wave obeys the following 
dispersion relation; 

^ _ k.v + ivk'^){ijj - k.v + ir^k^). (27) 

For a plane wave {uj real) and small damping rate, the 
wave vector describing the wave amplitude can be ap- 
proximated by: 



/c = fco + *"^^~~^0: 



(28) 



where fco 



VA+V 



is the wave vector in the absence of 



damping. The damping vector describing the energy or 
the wave action is then twice that of the wave amphtude: 



ftp — ; "-o- 



(29) 



In the stationary regime, the wave action is then damped 
as ([Jacques 19771): 



d_ 

dx 



[(v + VA)yl] = -fcdis [(v + va)-4] 



(30) 



This last equation will be used in Section 14.3.21 



4.3. Maximum amplification and entropy creation 

4.3.1. Simple estimate 

The maximum energy flux of the Alfven wave can 
be estimated by stating that amplification and dissipa- 
tion compensate exactly: |fcamp| = l^disl- Using Equa- 
tions (j25p and (j28l) , the wave number /cmax corresponding 
to the maximum amplitude can be estimated: 



h — 



1 



(31) 



where we implicitly assumed that the maximum amplifi- 
cation happens close to the Alfven surface and v + va ^ 
VA ~ VAo- In order to make a simple estimate, let us 
assume that the Alfven wave evolves in two steps: the 
energy flux is amplified without any dissipation in a first 
step, reaches a maximum i^max, and is instantaneously 
dissipated in a second step. Using Equation (|23l) . the 
maximum energy flux is then: 

-Fin VAin k'ni ^Ain k'mXd'u 



where the subscript "in" refers to the incident Alfven 
wave and should be replaced by out in the case where 
the Alfven wave is sent from below. The entropy created 
by the dissipation of the amplifled Alfven wave can be 
estimated with this maximum energy flux: 



1 



(33) 



Interestingly the maximum energy flux as well as the 
created entropy both diverge when the dissipation is 
decreased to zero, scaling like the square root of the 
Reynolds number. This suggests a singular behavior and 
large effects of the Alfven surface in a flow with low dis- 
sipation. The nonlinear saturation of the amplification 
at low dissipation is studied in Section 6. 




Fig. 5. — Amplification of an Alfven wave approaching tlie Alfven 
point: the transverse magnetic field ^/bJ^KB^' jg shown as func- 
tion of X. The black line is the result of the simulation (at the 
end of the simulation when a steady state is reached), the dotted 
line is the prediction by the ideal WKB analysis (Equation I I22II ). 
the dashed line is the WKB prediction that includes the effect of 
dissipation (deduced from Equation l|36|l ^. 

4.3.2. Exact calculation 

A more precise calculation can be done in the regime 
where all the dissipation takes place close to the Alfven 
surface, so that the velocity profile can be approximated 
by a linear function: 

V + VA — —Lu^x. (34) 

The damping vector can then be written explicitly as a 
function of x using Equation 



kdis — 



which allows us to solve explicitly Equation ([50)1 : 



(v + Va)^ = (Vin + VAin).4i, 



(35) 



(36) 



The corresponding transverse magnetic field profile is 
compared successfully with the simulations in Figure [5] 
(dashed line versus full line). In particular, the maximum 
magnetic field is: 



2 ^Ao 5 Bin 

6 ^Ain ^in^^dii 



(37) 



Using these explicit profiles, the amount of entropy 5S 
created can be computed with the ohmic dissipation as: 



5S^ 



+ 00 



■n 



HoPv 



(V X B)^ dx 



1 



n 

(38) 

which is very close to the simple estimate obtained in 
Equation ([55)1 . Note that in the case of a linearly polar- 
ized wave, the above equation should be multiplied by a 
f actor 1/2. 

IWilliami (|1975[ ) called the phenomenon described in 
this section an instability because the energy of the per- 
turbation (the incoming Alfven wave) increases exponen- 
tially. We prefer to call it an amplification, because it 
critically depends on the properties of the incident Alfven 
wave. For example the maximum amplitude reached is 



7 



proportional to its amplitude, and if the source of inci- 
dent Alfven wave is stopped the perturbations gradually 
disappear (unless the displacement is kept to a non zero 
value, see Section 5). Thus contrary to usual instabilities, 
the growth of perturbations at an Alfven surface needs a 
continuous incident wave and cannot grow to non linear 
amplitude from some small initial perturbations: the non 
linear saturation studied in Section 6 happens only if the 
incident amplitude is large enough for given dissipation 
coefficients. 

4.4. Pressure feedback 

In this subsection we investigate the origin of the pres- 
sure waves that are emitted up and down when the 
Alfven wave propagates in the gradients. Some insight 
about their physical origin can be gained by considering 
the energy budget of the flow. As shown in Section 4.1, 
the energy flux of an Alfven wave in the WKB regime 
increases when it approaches the Alfven surface, which 
implies a coupling with at least another wave in order 
to ensure energy conservation. As the WKB criterion 
precludes a linear coupling, the energy conservation re- 
quires a nonlinear coupling which causes a pressure feed- 
back. Figure |4] shows that, after a growth period, the 
pressure feedback levels off and remains at a constant 
value, suggesting a non linear coupling to a perturbation 
at zero frequency. By contrast a linear process would 
have retained the frequency of the incident Alfven wave 
uj. As will be shown in the following Section, the pressure 
feedback in the low frequency regime can display a time 
dependence with a frequency either 0, uj or 2uj. The ab- 
sence of non zero frequencies in the high frequency regime 
can be interpreted by the presence of many wavelengths 
in the gradients, which pressure feedbacks interfere de- 
structively. 

In order to estimate quantitatively the pressure feed- 
back, we consider the energy and mass budget of the 
region — < x < Lv- In the stationary state where 
the pressure feedback is constant, the net energy and 
mass fluxes should vanish: these two equations are suf- 
ficient to determine the two unknowns of the problem 
(the amplitudes of the acoustic signals propagating up 
and down). We further assume that the amplification is 
large, i.e. that the maximum energy fiux of the Alfven 
wave is much larger than the incident energy fiux. Ne- 
glecting the energy flux due to the incident Alfven wave, 
the energy and mass flux perturbations are due to: 

(i) the entropy wave, described by: 



follows: 



=1-7 — , 

P 7 

SPs = 0, 



(39) 

(40) 
(41) 



(ii) the acoustic wave propagating up (termed -) and 
down (termed -I-), described by: 



SP± 

pe- 



rn 

(43) 



The conservation of mass and energy can be written as 



(1 + '-)^ = (1 - ^)^^ + (1 - 7) -(44) 

Cin Pin ^ont Pout 



(45) 



Combining these two equations, the pressure perturba- 
tion is obtained as a function of the created entropy: 



(5P™p = 5P. 



SP„f ^SP,=- 



Cput -(7^1) Coutfout 
(Cin + ^^in) (cin + ^^^^^^ 

Cout + (7- l)Ci„Wi, 



(Cout - t'out) (cout + ^j^) 



PinSS, (46) 
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Using Equation (p8| . the pressure feedback can then be 
expressed as a function of the incident Alfven wave am- 
plitude and the other parameters of the problem. This 
analytical result is compared with the simulations in Fig- 
ure [6] The dependence on the amplitude and frequency 
of the incident wave, as well as the resistivity and the 
Mach number of the flow are reproduced very well. A 
slight discrepancy appears at very high resistivity be- 
cause the assumption of large amplification is not ful- 
filled. For incident waves with a large amplitude, the 
amplification saturates when the assumption of linear 
perturbation breaks down: this non linear regime is fur- 
ther studied in Section 6. In addition, a different regime 
appears at low frequency. This regime is studied in the 
following section. 

5. THE LOW FREQUENCY REGIME 

In this section we study the low frequency regime, 
where the WKB criterion is not satisfied: e = wv/w ^ 1. 
The fiow can be considered as stationary at each instant, 
because the period of the incident Alfven wave is much 
longer than the response time of the gradients. With this 
steady state assumption an analytical study of this fiow 
becomes possible in the vicinity of the Alfven surface. As 
previously, the transverse magnetic field is assumed to be 
arbitrarily small, and the velocity profile is described by 
the lin ear fu nction of Equation p4|) . As shown in ap- 
pendix IB.H the transverse magnetic field then follows 
the second order ordinary differential equation: 



2 



ox 



y 



6B, 



(48) 



Note that the z component follows the same equation. 
The solution that does not diverge for large a; is a Gaus- 
sian function with a typical width a;dis: 



SBy = 6B„ 



(49) 



This suggests that in this regime part of the displacement 
of the field lines could accumulate in a region of width 
Xdis around the Alfven surface. This is confirmed by the 
simulations, as shown in Figure [T] The displacement of 
the field lines induced by this profile can be computed 



as: 
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Fig. 6. — Entropy creation and pressure feedback: comparison between the simulations and the analytical calculation. The four plots 
show the dependence on the following parameters (from top to bottom and left to right): the amplitude of the incident Alfven wave SB^^, 
its frequency u, the resistivity represented here by the parameter x^^c^ = ^ v/'^V i s-i^d the sound speed cq. As each of these parameters is 
varied, the others are kept constant to their baseline values. The entropy creation 5S (normalized by (Bo/^in)^ in order to be independent 
of the incident wave amplitude in the linear regime) is shown with crosses for the simulation result, and dashed lines for the analytical 
prediction. The upstream pressure feedback <5Psup, normalized by the magnetic pressure of the incident wave SBf^, is represented by 
triangles (simulation results) and full lines (analytical prediction). Finally, in the bottom right plot, the pressure feedback downstream is 
shown by crosses (simulations) and dotted line (analytical) with the same normalization. Notice the two frequency regimes separated by a 
vertical dotted line at ujuj^ = 1/-k. The agreement between analytics and numerics is excellent, except at large amplitude of the incident 
wave where non linear effects cause a saturation of the pressure feedback and entropy creation (see Section |6] for more details). 

Figure[7]compares the corresponding profile (dashed line) 
with the simulations (full lines). This confirnis the Gaus- 
sian shape and the frequency dependence (contained in 
the normalization) in the low frequency limit. Further- 
more Figure [7] shows that a low frequency Alfven wave 
behaves like a stationary displacement of the field lines 
(here induced by sending only half a wavelength of a 
linearly polarized Alfven wave). As shown in Figure [7] 
however, Equation ([52]) overestimates (5i?max by ^ 10%. 
This can be interpreted by the fact that part of the dis- 
placement of the incident wave goes into a downward 
propagating Alfven wave created by linear coupling in 
the gradie nts. A more elaborate estimate described in 
Appendix IB. 21 resolves this slight discrepancy. The en- 
tropy created by this circularly polarized wave is constant 
and estimated with Equations ([5T|) and (|52p . which can 
be written in a similar form as Equation ([38]): 



The entropy created by this profile is: 



The relation between the incident Alfven wave ampli- 
tude and ^Bmax can then be determined approximately 
by assuming that the downstream field lines are not per- 
turbed at all. This amounts to neglecting the downward 
propagating Alfven wave created by linear coupling of 
the incident Alfven wave in the gradient, and leads to 
an error of ~ 10%. This simple approximation allows 
to match the transverse displacement of the incident low 
frequency Alfven wave with that contained in the vicin- 
ity of the Alfven surface (Equation ([501 ) and deduce the 
value of (5i3inax- 

In order to check the above description, let us first 
consider a circularly polarized Alfven wave which dis- 
placement is symmetric with respect to the initial flow. 
The displacement of such a wave is constant and equals 
^/&\p^r^ = SBin/ikinBo), giving: 



6S 



1 Fn 



2ttPqVo 



(53) 



SB. 



(52) 



This formula is compared successfully with the simula- 
tions in Figure [6| (left bottom panel). 
In the case of other Alfven wave polarizations, the dis- 
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Fig. 7. — Left panel: Magnetic field fines for a circularly polarized Alfven wave in the low frequency regime (e = u)^/uj = 15). Right 
panel: Transverse magnetic field profile (normalized by the displacement amplitude of the incident wave) for different frequency Alfven 
waves in the low frequency regime compared to the profile produced by a stationary displacement of the field lines. From bottom to top 
(and black to red in the online version): a stationary displacement, e = 15, e = 6, e = 3, and the three higher cur ves (e = 1.5, 0.6 and 0.3) 
which arc not in the low frequency regime. The dashed line shows the analytic estimate given by Equation l|49|l and ((52j, which has the 
correct shape but overestimates the maximum magnetic field by ^ 10%. 



placement of the field lines (and therefore the entropy 
created) can be time dependent. However the low fre- 
quency assumption ensures that this timescale is much 
longer than the dynamical timescale, such that the os- 
cillation can be considered as a succession of station- 
ary states and the analysis of Section 4.4 applies. The 
pressure feedback is proportional to the entropy creation 
(Equations P5|) - (|T7|) ) and therefore to the square of the 
displacement. As a consequence the time dependence 
of the pressure feedback can be understood from the 
time dependence of the displacement. The trajectory 
of a field line and the corresponding pressure feedback 
is illustrated in Figure [8] in four different cases. For a 
circularly polarized Alfven wave started as described in 
Section 2.2, the field line describes a circle around its 
initial position and thus creates a stationary feedback 
(i.e. at zero frequency). If the first quarter of a period is 
omitted (dashed line), it describes an off-centered circle, 
leading to a displacement and a pressure feedback that 
oscillates with the same frequency as the incident wave. 
Similarly, linearly polarized Alfven waves show two dif- 
ferent types of time dependence: the centered trajectory 
leads to an oscillation of the displacement in the range 
~5ya < 6y < Syo and a pressure feedback with a doubled 
frequency. The off-centered trajectory oscillates between 
< Sy < 26yo, thus the pressure has the same frequency 
as the incident wave. The maximum displacement of off- 
centered waves is twice that of centered waves, leading 
to a maximum pressure feedback that is approximately 
four times larger (Figure [5]). 



Magnetic field lines trajectory 




linear polarization : 

centered 

- - - off-centered 

circular polarization : 

centered 

off-centered 




6. 



SATURATION OF THE AMPLIFICATION AT LARGE 
AMPLITUDE OR LOW DISSIPATION 



When the amplitude of the incident Alfven wave is 
increased or when the dissipation is decreased, the maxi- 
mum amplitude of the Alfven wave increases and eventu- 
ally becomes non linear, and the analysis of the preceding 
sections no longer applies. This non linear regime starts 
when the Alfven wave significantly changes the station- 
ary flow profiles (e.g. of pressure, density, velocity), i.e. 
when it becomes compressible. In this section we de- 
scribe this non linear dynamics and show how it affects 



time/TjjK 

Fig. 8. — Time dependence of the pressure feedback in the simu- 
lations in the low frequency regime (here a; = 5 X 10"^, e = 3). The 
black curves correspond to circularly polarized Alfven waves, the 
thick gray (red in the online version) curves to linearly polarized 
waves. The displacement is either symmetric with respect to the 
initial flow (full lines), or one sided (dashed lines). 



the Alfven wave amplification and the creation of the 
pressure feedback. To explain the numerical results, we 
give an approximate analytical description of the satu- 
ration, which takes into account the non linearity of the 
Alfven wave but not of the pressure feedback. The still 
higher amplitude regime where the pressure feedback has 
a significant impact on the flow cannot be studied with 
our numerical set up, because the acoustic wave interacts 
with the grid boundaries in a non physical way. 

The dynamics of an Alfven wave in the linear regime is 
very similar in the different configurations: whether it is 
sent from above or from below gives the same amplifica- 
tion and pressure feedback, the only difference between 



10 



0.5 
0.0 
-0.5 r 

-1.0 : 



- linear pol. 
circular pol. 





Fig. 9. — Magnetic field along the y direction (up) and prop- 
agation velocity of the Alfven wave vx+vA (bottom) at the end 
of the simulation. The Alfven waves are sent from above with an 
amplitude of (5 Bin = O.ISBq. The full line corresponds to a lin- 
ear polarization, the dotted line to a circular polarization, and the 
dashed line to the unperturbed stationary flow. The linearly po- 
larized wave steepens and therefore dissipates efficiently. On the 
other hand the circularly polarized wave does not steepen, reaches 
larger amplitudes and pushes the Alfven surface downward. 



the circular and linear polarization is the time depen- 
dence in the low frequency regime. On the contrary, the 
non linear dynamics described in this section differs for 
these configurations, we thus describe them separately in 
the next subsections. 

6.1. Steepening of linearly polarized Alfven waves 
The propagation velocity of an Alfven wave non lin- 



early changes as described hy v + va in Equation (|A6[) 
Thus the propagation velocity of a maximum of trans- 
verse magnetic field pressure is different from that of a 
minimum, leading to a steepening of the wave (Figure |9]). 
The amplitude above which this steepening happens can 
be estimated with the criterion that the minimum and 
the maximum (separated by one quarter of a wavelength) 
have time to catch up before the wave reaches the Alfven 
surface: 



5'^(v + va)- 



VA 



n V + VA 
2 uj ■ 



(54) 



Using Equations ([^ and one can get the corre- 

sponding transverse magnetic field strength for a given 
incident amplitude: 



6B„ 



„ a^v(«in +fAin) ,2 2 N.p 

W^WAOWAin 



1/3 



(55) 

The steepening of the wave damps it efficiently, thus 
Equation (I55|) is a good estimate of the maximum mag- 
netic field reached by a linearly polarized Alfven wave in 
the non linear regime (Figure 1101 upper panel) . As in 
Section 4.3.1, we then use the corresponding maximum 
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Fig. 10. — Saturation at large amplitude of linearly versus circu- 
larly polarized Alfven waves sent from above. The maximum trans- 
verse magnetic field (upper panel) and the entropy created (lower 
panel) are plotted as a function of the incident wave amplitude. 
Both quantities are normalized to be constant in the linear regime, 
so that the onset of the nonlinear regime is clearly visible. Symbols 
are the simulation results (diamonds for linear polarization, plus 
signs for circular polarization), while lines represent the analyti- 
cal estimates (full (dotted) lines for linear (circular) polarization). 
The two line segments represent the linear solution (horizontal) and 
nonlinear estimates (note that for the circular polarization this is 
only an estimate of the onset of the nonlinear regime, unlike the 
linear polarization). 

flux to estimate the amount of entropy created: 

^ max 
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fJ-aPoiCo - fAo)("m + VAin) SB^^^ (57) 



(56) 
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This estimate gives the good scaling with SBm but is a 
factor 2 smaller than the result of the simulations (Fig- 
ure [ini lower panel) . This is reasonable given the sim- 
plicity of the argument. 

6.2. Circularly polarized Alfven waves 

Circularly polarized Alfven waves do not steepen like 
the linearly polarized ones, because their magnetic pres- 
sure does not oscillate over a wavelength. Even at non 
linear amplitude, they propagate without deformation 
in a uniform medium since they are an exact solution 
of MHD equations. In a non uniform medium however, 
this is no longer true and the amplification of the Alfven 
wave creates a gradient of v + va that is non linear (using 
Equations dHJ and (j25|)): 



V {6'^ {v + va)) ^ 



VA 



SB^ 



V + va2p (c^ 



(58) 



This effect is significant if the gradient created by the 
wave is comparable to the stationary gradient: 



(59) 



Using Equations ((58)) and ([22]) for a given incident 
amplitude, the corresponding transverse magnetic field 
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Fig. 11. — Same as Figure |9] but for a circular Alfven wave sent 
from below with an amplitude SBout/Bo = 0.18 (dotted lines), and 
a stationary displacement of the field lines (full line) which is repre- 
sentative of the low frequency regime. The stationary displacement 
is induced by sending half a wavelength of a circularly polarized 
Alfven wave of frequency uj = 1 and amplitude 5Bin/Bo = 0.5 
(the asymmetric wave of Figure [Sj , the result depends only on its 
displacement 2SB[^/ (Boki^). 



strength is obtained: 
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(60) 



Comparing this with the maximum transverse magnetic 
field in the hnear regime (Equation (|37p ). gives the 
threshold of the non linear regime. The behavior is dif- 
ferent if the wave is sent from above or from below. The 
gradient created by an Alfven wave sent from above is 
opposed to the stationary gradient (because v + va < 
at X > 0). Therefore it weakens the gradients and pushes 
the Alfven surface downward, as illustrated in Figure [9l 
The weaker gradients lead to a less efficient amplification 
and pressure feedback (Figures [51 and [TU)) . 

In contrast, the gradient induced by a wave sent from 
below steepens the stationary gradient (because v + va > 
at X < 0). This causes more amplification which causes 
more steepening, leading to the build up of a discontinu- 
ity (Figure ITlT) . Across this discontinuity the flow is ac- 
celerated to a velocity close to (though somewhat slower 
than) the Alfven speed, while the Alfven wave is ampli- 
fied. Just above the threshold, the discontinuity is weak 
and close to the Alfven surface. As the amplitude is in- 
creased, the discontinuity becomes stronger and moves 
away from the Alfven surface. Contrary to the case of 
an Alfven wave sent from above this nonlinear effect does 
not diminish the amplification, thus no saturation occurs 
for SBin/Bo < 0.3 (Figure [T2|) . For larger amplitudes the 
discontinuity becomes so strong that it reaches the lower 
grid boundary, leading to non physical effects. 

6.3. The low frequency regime 

As described in Section 5, a low frequency wave in 
the linear regime accumulates in a region with a width 
Xdis centered on the Alfven surface. In the non linear 
regime, the magnetic pressure accumulated around the 
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Fig. 12. — Same as Figure lTOl but for a circular Alfven wave sent 
from below (plus signs, dotted lines) and a stationary displacement 
representati ve o f the low frequency regime (diamonds, full lines). 
As in Figure [TT1 the stationary displacement is induced by sending 
half a wavelength of a circularly polarized Alfven wave of frequency 
ui = 1 and amplitude SBi^ (the asymmetric wave of Figure [SJ , the 
result depends only on its displacement 2S Bi^ / {Boki^) ■ 

Alfven surface flattens the velocity proflle. Eventually 
the Alfven wave accumulates over a wider region veri- 
fying V + VA — and located downstream of the Alfven 
surface (Figure [TT]). The transverse magnetic field proflle 
can be estimated with the condition v + va = 0: 

-ujvx + S'^{v + va)-0, (61) 

thus using Equation (|A6|) : 
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(62) 



This proflle is compared with a simulation in Figure 1111 
The (relatively small) discrepancy can be attributed to 
the effect of the pressure feedback on the velocity proflle 
which was neglected. The width xinf of the region can be 
approximated as in Section 5 by matching its displace- 
ment with that of the incident Alfven wave: 
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thus giving the maximum magnetic field: 
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(63) 



(64) 



The entropy created by this proflle can be roughly esti- 
mated using: 



SS = 



poPv 



(V X B)^ dx 



r]SB'^ 



(65) 



where L is the scale over which SB varies. This entropy 
creation is dominated by the sharp decrease of SB located 
near —Xmi- The length scale L can be roughly estimated 
assuming that the width of this sharp decrease is due to 
the dissipative scale. The width is then such that this 
diffusion speed due to the resistivity (typically v^is 
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Fig. 13. — Saturation at low dissipation in the low frequency 
regime (represented by a stationary displacement, i.e. half a wave- 
length of a = 1 wave) . The maximum transverse magnetic field is 
plotted as a function of the dissipative scale xjis, for three values 
of the incident wave amplitude: SB\^/Bq = 8.10"'* (plus signs), 
(5Sin/_BQ = 8.10-3 (diamonds) and SBin/Bo = 3.1Q~^ (triangles). 
The dashed line represents the analytical estimate in the linear 
regime, while the two full lines correspond to the nonlinear esti- 
mate for SBin/Bo = 8.10-3 (up) and SBi^/Bq = 3.10-^ (down). 

■^) compensates the velocity at which the Alfven wave 
propagates: Vdis ~ i-^vXinf- Therefore, we get: 



L 



(66) 



Finally, using the four previous equations, the entropy 
created can be expressed as: 



Po(c§- wio) / 9 
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(67) 



Note that a similar argument in the linear regime gives 
the same scaling as the more precise calculation of Sec- 
tion 5. The estimates of the maximum magnetic field 
and entropy creation are compared successfully with the 
simulations in Figure 1121 

While the linear amplification diverges at low dissipa- 
tion, it is important to mention that Equations ()64|) and 
(j67p show that the nonlinear saturated state is indepen- 
dent of the dissipation coefficient. This is checked with 
the simulations in Figure 1131 The same phenomenon was 
found in Section 6.1 for the saturation of linearly polar- 
ized waves in the WKB regime (Equations (|55|) and (f57| ). 
This can be interpreted by the fact that the amplitude 
at which the Alfven wave is dissipated is controlled by 
nonlinear effects rather than by the dissipative scale. 

7. CONSEQUENCES FOR CORE COLLAPSE SUPERNOVAS 

We now consider the consequences of an Alfven sur- 
face in the context of core collapse supernovae, focusing 
on the phase of stalled accretion shock when matter from 
the outer layers of the core accretes on the forming proto- 
neutron star. What maximum amplitude do the Alfven 
waves reach? Can the pressure feedback significantly af- 
fect the shock dynamics? Using the analytic formula 
derived in this paper, we estimate the conditions under 
which these effects could be significant. 

7.1. Origin of the Alfven waves: 

Several instabilities cause non radial motions, which 
can potentially be a source of Alfven waves: both neu- 
trino driven convection and the standing accretion shock 
instability (SASI) operate between the neutrinosphere 
and the shock, while Ledoux convection takes place below 
the neutrinosphere. Powerful Alfven waves could also be 



created by the g-modes of the PN S if these are ex c ited t o 
large amplitudes as sugg ested bv iBurrows et all (|2006D . 
Below we estimate the typical amplitude and frequency 
of the Alfven waves created by these three phenomena. 

SASI in a magnetized flow was studied linearly by 
IGuilet fc Foglizzol (|2010D . Although SASI can be sta- 
bilized by a radial fleld, this stabilization effect is ex- 
pected to be less acute if the field is oblique, since the 
horizontal component ca n enhance the growth of SASI. 
IGuilet fc Foglizzol ([2010) showed that the vorticity cre- 
ated by the shock oscillations is distributed into slow and 
Alfven waves. The proportion of these two types of waves 
depends on the field geometry and the wave vector of the 
SASI mode. Except for very particular cases (for exam- 
ple in a purely radial magnetic field, only slow waves 
are created) one can expect that a significant fraction of 
the vorticity is in the form of Alfven waves. Further- 
more, downward and upward propagating Alfven waves 
are created with roughly the same amplitude (strictly 
in the weak field limit). As a first estimate we will 
thus assume that the amplitude of the upward propa- 
gating Alfven wave of interest for the Alfven surface is a 
few times smaller than the maximum transverse veloci- 
ties induced by SASI. According to published numerical 
simulations, the deformations of the shock induce trans- 
verse motions with a sp eed comparable to or even greater 
than the sound speed fBlon din et al.|[2"003t iScheck et al.l 
2008; Marck & Janka 2009). It is thus reasonable to take 
an amplitude of 5v ^ lO'^km.s"^, which corresponds to 
Sv/c ^ 0.3. The osciUation period of SASI is deter- 
mined by the advection time from the shock to a radius 
of strong deceleration slightly above the neutrinosphere 
(jFoglizzo et alJ 120071: IScheck et ai]|2008D . and is usually 
in the range: TAif ~ 30 - 50ms (w ~ 100 - 200 s^^). 

The production of Alfven waves by Ledoux con- 
vection in t h e pro to-neutron star was considered by 
iSuzuki et al.l ()2008D . The wave amplitude may be es- 
timated by the typic al velocities in the PNS convection : 
fa ~ 3 X lO '^km.s-i (|Keil et al.lll996t IBuras et al.ll2006l: 
IDessart eFa l. 200(3,). The period may be roughly esti- 
mated with the overturn time of a convection cell with a 
si ze I ^ 10km: TAif ^ l /v ^ 3ms (w - 2.103 s-^). 

IBurrows et al.l (120061 ) witnessed powerful g-mode os- 
cillations that had a typical period of T ~ 3 ms and 
a displacement amplitude of Ar ~ 3 km. The Alfven 
waves created by these oscillations may be estimated to 
have similar amplitude and frequency : Sv ~ uiAr ~ 
6000 km. s^^. However one should note that such large 
amplitudes are still controversial since they have not 
yet b een confirmed by other groups (iMarek fc Janlaj 
120091). Furthermore iWeinberg fc QuataertI H20Q8f ) sug- 
gested that non linear coupling with short wavelength 
modes unresolved in the simulations could saturate this 
g-mode at a much lower amplitude. 

These amplitudes in the range Sv ^ 0.3— 1.10'* km. s"^ 
are a significant fraction of the sound speed c ~ 3.10** km 
{6v/c ~ 0.1 — 0.3). Therefore one would expect such 
waves to be in the nonlinear regime studied in Section 6. 



7.2. Amplification time scale: 

The timescale of the Alfven wave amplification and 
pressure feedback is controlled by the stiffness of the ve- 
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locity gradient: 



V 



(68) 



where iv is the typical size of the gradient and v a typ- 
ical advection speed. This timescale depends crucially 
on the radius of the Alfven surface. If the Alfven sur- 
face is inside the proto-neutron star, the advection speed 
due to the contraction of the PNS is very slow and the 
timescale could become too long to affect the dynamics 
of the supernova explosion. 

The situation would be more favorable if the Alfven 
surface lies outside the PNS. For example, iScheck et al.l 
(|20081 showed that the velocity gradient peaks at a ra- 
dius slightly above the neutrinosphere and reaches val- 
ues of 100 — 300 s~^ (their Figure 14). This translates 
into a quite fast growth timescale: r ^ 3 — 10 ms. This 
timescale, associated to the advection time across the 
deceleration layer with a size Lv ~ 10 — 20 km, is sig- 
nificantly shorter than the period of SASI, which corre- 
sponds to the advection time from the shock to the PNS 
surface on a distance of Tsh — T-pArs ^ 100 km. As a conse- 
quence, SASI-created Alfven waves would lie marginally 
in the low frequency regime (wv/w ~ 1 — 3), while those 
originating from the PNS convection would lie in the 
WKB regime (wy/w ~ 0.05 - 0.15). 

We conclude that the Alfven wave amplification can 
be faster than other dynamical timescales, if the Alfven 
surface lies in the strong gradients slightly above the neu- 
trinosphere. In the next subsection, we discuss whether 
this corresponds to a realistic magnetic field strength. 

7.3. Magnetic field strength and radius of the Alfven 
surface: 

Because the infall velocity decreases to zero at the cen- 
ter of the proto-neutron star, an Alfven surface should 
be present in the flow regardless of the magnetic field 
strength. However, its consequences and their signifi- 
cance for core collapse supernovae depends on the Alfven 
surface radius. This radius is controlled by the magnetic 
field strength, which is unfortunately extremely uncer- 
tain. A very weak magnetic field would place the Alfven 
surface close to the center inside the proto-neutron star, 
while stronger field would move it outward to larger radii 
where the infall velocity is larger. 

As discussed in the previous subsection, if the Alfven 
surface is inside the PNS the amplification timescale may 
be too long to affect the dynamics before an explosion 
sets in. Therefore we concentrate on the possibility that 
the Alfven surface lies above the PNS surface. What 
strength of the magnetic field is required for this to hap- 
pen? Assuming an advection speed of u ^ 1000 km. s~^ 
and a density oi p ^ 10^° g.cm^'^, yields a magnetic field 
at the Alfven surface of i? ~ 3.5 x 10^'^ G. Compared 
with other supernovae studies, this corresponds to an in- 
termediate strength. Indeed, magnetic models of super- 
novae explosions often find or assume a much stronger 
field of ~ 10^^ — 10^^ G with a m agnetic pressure com- 
parable to the thernial pressure (lAkivama et al.l [20031: 
Shibata et all 120061: [Burrows et al .l 120071: iSu zuki et al.l 
2008t lObergaulinger et al.l I2009D . while here the mag- 
netic pressure remains much smaller than the thermal 
one: Pma,g/P ~ 7/2 x v\/c^ < 10"^. On the other hand. 



this magnetic field strength is significantly larger than 
th e extrapolat i on fro m the stellar evolution calculation 
of lHeeer et all ((20051) . 

In order to get an idea of the plausibility of such a field, 
it is instructive to estimate the magnetic field of the re- 
sulting neutron star and to compare it with observed val- 
ues of pulsars and magnetars. Assuming that the mag- 
netic flux is conserved (with either the usual prescription 
Bp^/^, or Br'^ for a radial fleld) and extrapolating the 
magnetic field at the Alfven surface to the final density 
and radius of a cool neutron star {p ^ 3.10^* g.cm^^, 
r ~ 10 km) yields B - lO^^ - lO^^ G. While this is much 
larger than usual pulsars, it corresponds to the range of 
magnetic field strength typical of magnetars. Therefore, 
we conclude that the Alfven surface may play an im- 
portant role in supernovae leaving a magnetar behind, 
but should not be significant in the more common su- 
pernovae forming normal pulsars. One should however 
keep in mind that the correspondence between the mag- 
netic field of the PNS and that of much older neutron 
stars is very uncertain, because the magnetic field might 
decrease during the super novae explosion and la ter evo- 
lution, as hypothesized by I Wheeler et al.l (j2002D . If this 
were the case, magnetic effects would play a role in a 
larger subset of supernovae. We should also caution that 
our description of the Alfven surface dynamics applies 
only if the m agnetic field comes from the progenitor star 
(IFerrario fc Wickramasin g hc 2006.) , and not if created in 
situ by a dy namo or the tapping of the r otational energy 
o f the PNS (iThompson fc Duncanl ll99l. 

ISuzuki et al.l ( 2008[ ) also studied the effect of Alfven 
waves on the dynamics of core collapse supernovae. They 
assumed a magnetic field of i? ^ 1 — 3.10^^ G at the 
surface of the proto-neutron star (r = 50 km) , which is 
one to two orders of magnitude larger than considered in 
this paper. As a consequence, their Alfven surface lies in 
the supersonic flow above the shock, while we restricted 
our study to the case where the Alfven surface is in the 
subsonic part of the flow. 

7.4. Viscosity: 

In the linear regime, the maximum amplification and 
the pressure feedback critically depend on the dissipation 
scale (due to viscosity or resistivity). However, the large 
amplitude of the incident Alfven wave would probably 
place it in the non linear regime discussed in Section 6, 
in which the dynamics is independent of the dissipation 
coefficient. Therefore the dissipation should not play an 
important role unless it is much larger than considered 
in this paper ; we check below that this is not the case. 

The hot nuclear fluid constituting the proto-neutron 
star is an extremely good conductor but is rather vis- 
cous due to the diffusion of neutrinos (resulting in a 
huge magnetic Prand tl number Pm ~ i^/jy ^ 10^^, e.g. 



The dominant diffusion 



IThompson fc DuncanI 
therefore comes from the neutrino viscosity, which can be 
estimated in the region optica lly thick to neutrinos u s- 
ing the diffusion limit (see e.g. IThompson et al.l ()2005D ): 
1/ ~ lO^^cm^.s— 1, slightly below the neutrinosphere 
where it is highest. This translates into a dissipation 
scale Xdis ^ \/ J^/wv ~ 0.2 km and a typical Reynolds 
number Re ^ (Ly/xdis)'^ ~ 10** (using wv ^ 200 s~^, 
Lv ~ 20km as suggested in the previous paragraph). 
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This value is in the range explored in this paper and 
yields a significant amplification of an Alfven wave be- 
fore it is dissipated. 

The steep gradients proposed to be a good place for a 
strong effect of the Alfven surface lie above the neutri- 
nosphere. There, the matter is optically thin to neutri- 
nos and this neutrino viscosity is unfortunately not well 
defined. However, interactions with streaming neutrinos 
sti ll have a similar d amping effect, which was estimated 
by lThompson et aD ((2005) by a damping rate F ^ 30 s~^. 
As this is smaller than the growth rate of the amplifica- 
tion, we do not expect the neutrinos to be able to prevent 
the amplification very efficiently. 

7.5. Amplitude of the pressure feedback: 

In order to get a simpler scaling of the pressure feed- 
back, we assume that at the Alfven surface, the fiow is 
very subsonic, that is: v ^ va c. Equation (j46|) giving 
the pressure feedback propagating toward the shock then 
simplifies into: 



SP 



— - -^5S > SS. 



(69) 



Since Alfven waves created by convection and g-modes 
in the proto-neutron star lie in the WKB regime, we use 
the formulation of Section 6.1. Under the same assump- 
tions, Equation (|57)) simplifies into: 



SP_ 



7 (2t: 



/ «in+t^Ain faV^^ 



V^AinWo C 



(70) 



This gives 5P/P ~ 0.05 for PNS convection (using 
the values: uj^jio ~ 0.1, bvjc ~ 0.1 and (win -I- 
t'Ain)/V^Ainfo ^ 1) and 6P/P ^ 0.1 for g-modes (us- 
ing (5w/c~ 0.2). 

The Alfven waves created by SASI lie marginally in 
the low frequency regime, therefore we use the nonlin- 
ear formulation developed in Section 6.3. Equation (|67|) 
simplifies into: 

The values u)\j /uj ^2, (v -\- va)/va ^ 1 and 5v/c ^ 1/3 
give a very significant pressure feedback: 6P/P ^ 2. The 
corresponding maximum transverse magnetic field can be 
estimated with Equation ([Ml) : 



SBn 



^UJyV + VA 

DC OV 

LO VA 



1/3 



(72) 



yielding (5i?max ^ 1-5 x 10^^ G. Of course this very large 
value of the pressure feedback should not be taken at face 
value, since our analysis assumes that the pressure feed- 
back does not significantly modify the upstream flow. A 
self consistent calculation should consider the shock as 
the source of Alfven waves in order to capture the effect 
of the pressure feedback on the Alfven wave injection. 
We also caution that such a large pressure feedback cor- 
responds to maximum values of the transverse magnetic 
field with a pressure comparable to the thermal pressure, 
i.e. much greater than the longitudinal field. Whether 
this can really happen should be checked further in par- 
ticular with multidimensional simulations, which might 



show a different dynamics that could saturate the ampli- 
fication earlier that in ID. 

We conclude that if the Alfven surface lies in the strong 
gradients above the neutrinosphere, the amplification of 
the Alfven waves created by SASI generates a significant 
pressure feedback. SASI is more efficient than the PNS 
convection and g-modes because its Alfven waves have a 
larger amplitude and a lower frequency. A linear acoustic 
feedback is known to exist in the advective-acoustic cy- 
cle, which is believed to be the mechanism driving SASI 
(jFo glizzo e t al. 2007; Schcck ct al. 2008). By contrast, 
this new pressure feedback comes from a non linear cou- 
pling. As a consequence, it is potentially important only 
in the non linear phase of SASI. Also the frequency is 
not necessarily the same as that of SASI, but has Four- 
rier components at the frequencies 0, Wgasi and 2ti;sasi- 
More importantly, instead of oscillating between nega- 
tive and positive values 5P is here always positive. It is 
thus expected to be always pushing the shock outwards. 

7.6. Energy budget 

For the processes described in this paper to have a sig- 
nificant impact on the dynamics, the Alfven waves should 
have access to an appropriate source of energy. As argued 
in Section 7.1, a fraction of the kinetic energy of SASI and 
PNS convection could be tapped. If this were the only 
source of energy, these Alfven waves would probably re- 
main subdominant. However, as shown in Section 4, the 
energy flux of Alfven waves increases as it approaches the 
Alfven surface. During this amplification they thus ex- 
tract an additional energy from the stationary fiow. The 
stationary flow is potentially a much greater energy reser- 
voir than the kinetic energy, if the thermal or the gravi- 
tational energy can be tapped. An important question is 
thus whether the Alfven waves can extract a significant 
fraction of the thermal energy or if it is limited to the 
kinetic energy. This depends on the nonlinear saturation 
of their amplification. In Section 6 we found that the am- 
plification stops roughly when the magnetic pressure is 
comparable to the thermal pressure, which suggests that 
the energy of the Alfven waves can become comparable 
to the thermal energy. One should note however that 
our simulations were limited to a regime of parameters 
where the kinetic and thermal energies are comparable. 
Additional simulations are needed to address this point. 

8. CONCLUSION 

The main findings of this paper may be summarized 
as follows: 

• A small amplitude Alfven wave propagating toward 
an Alfven surface is amplified until its wavelength is 
as small as the dissipative scale Xdis- The timescale 
for this amplification depends on the stiffness of the 
gradient: — d(v + VA)/dx. 

• The amplification of the Alfven wave creates a pres- 
sure feedback that increases the upstream pressure 
through a non linear coupling. This pressure feed- 
back is stationary for a high frequency incident 
wave, while for a low frequency wave it can also 
display a time dependence with either the same 
frequency or twice the frequency of the incident 
wave. 
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• The subsequent dissipation of the amplified Alfven 
wave leads to an entropy creation, which can be 
related to the efficiency of the pressure feedback: 
6P/P ^ 5S (in a very subsonic flow). 

• In the linear regime, the smaller the viscos- 
ity/resistivity, the more efficient the amplification 
and the consequent coupling because the dissipa- 
tion happens later on (i.e. at smaller scale). The 
relevant quantities (maximum amplitude, pressure 
feedback, entropy created etc..) scale as (77 + 
i/)~^/^, and could thus be very large in a weakly 
dissipative flow. 

• At large amplitude of the incident wave or at low 
dissipation, the amplified Alfven wave becomes non 
linear and affects the background flow. A linearly 
polarized wave then steepens, which leads to its ef- 
ficient dissipation and therefore to the saturation of 
the amplification. In the low frequency regime, the 
nonlinear saturation is mediated by the flattening 
of the velocity profile in the vicinity of the Alfven 
surface. 

• In this nonlinear saturated regime, the dynamics 
becomes independent of the dissipation coefficients, 
because the amplitude at which the Alfven wave is 
dissipated is controlled by nonlinear effects rather 
than by the dissipative scale. 

An application to core collapse supernovae leads to the 
following conclusions: 

• Powerful Alfven waves can be created by the proto- 
neutron star convection, the g-mode oscillations 
and the SASI induced shock oscillations. 

• The amplification of these waves near the Alfven 
surface can be fast enough to affect the dynam- 
ics, if the magnetic field is strong enough to place 
this surface in the region of strong deceleration just 
above the neutrinosphere. Assuming magnetic flux 
conservation, this would correspond to the forma- 
tion of a magnetar. 

• In this case, the pressure feedback is dominated 
by the SASI created Alfven waves, and could 
contribute significantly to the pressure below the 
shock. The magnetic field could be amplified to 
a strength of ^ 10^^ G. The additional pressure 
can be expected to push the shock to larger ra- 
dius, which might be favorable for the success of a 
potential neutrino driven explosion. 

Our setup is voluntarily chosen as simplified as possible 
in order to show the physical processes at work in their 
simplest form and to render analytical work tractable. 
Because of its simplicity, our analysis suffers from the 
following major limitations: 

• We restricted our study to a one dimensional ge- 
ometry. It is possible that multidimensional phe- 
nomena neglected here stop the amplification of 
Alfven waves at a somewhat lower amplitude than 
in one dimension. Multidimensional processes that 



coul d dissipate Alfven waves include phase mix- 
ing (|Hevvaerts &: Priestl I1983D . nonlinear multi- 
wave interactions which lead to a turbulent cas- 
cade ( iGoldreich & : Sridhar 199^, and parametric 
decay instability (Goldstcini il978l iDel Zanna et al.l 
[2001) . This last process can in principle also hap- 
pen in one dimension, although we did not see any 
evidence for it in our simulations. 

• The magnetic field was assumed to be parallel to 
the advection velocity, which is a very singular sit- 
uation (for example, no Alfven discontinuity, even 
non evolutionary, is possible under these circum- 
stances). This might be expected in a very sub- 
Alfvenic fiow where the field lines are able to force 
the gas to follow them, however there is no rea- 
son why this should also be the case in the super- 
Alfvenic flow. 

• Incident Alfven waves were imposed at the bound- 
ary condition, without specifying how they were 
created. Strictly speaking, SASI cannot be invoked 
in the oversimplified field geometry assumed, be- 
cause SASI would create only slow waves but no 
Alfven waves. Alfven waves would be generated by 
SASI only if the field were oblique. 

• We studied only the behavior of pure sinusoidal 
waves, which are meant to model SASI induced os- 
cillations. This approach is a crude approximation 
since the fiow is known to be turbulent. The inter- 
action of such a flow with the Alfven surface should 
be the subject of future multidimensional simula- 
tions. 

• We used a simple adiabatic equation of state. Ne- 
glecting the cooling by neutrinos renders our sta- 
tionary pressure profile unrealistic. 

• We assumed a planar geometry. While this is jus- 
tified in the vicinity of the Alfven surface (where 
most of the dynamics takes place) , the convergence 
should be taken into account for a self consistent 
description from the source of Alfven waves to their 
amplification. 

• We assumed a large scale magnetic field that 
connects the region of emission and the Alfven 
surface. These two regions may be mag- 
netically disconnected if the field is turbulent 
(iThompson fc Duncan' Jl993l : lObergaulinger et al.l 
|2009; Endevc et al. 201C 

• Last but not least, we took as a free parameter the 
strength of the background magnetic field without 
explaining its origin. We thus implicitly assumed 
the field to originate from the progenitor star, and 
our results directly apply only in this scenario. The 
magnetic field already present in the collapsing core 
might or might not be sufficiently strong to place 
the Alfven surface above the neutrinosphere. If it 
is not strong en ough, it could still be amp lified in 
situ by the MRI (lObergauhnger et aLll2009 ') or by a 
dynamo process (|Thompson fc Duncan| [l993). but 
the extension of such a field beyond the Alfven sur- 
face is uncertain. 
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A realistic description of the Alfven surface dynamics 
will therefore require further investigation, in particular 
through multidimensional simulations. One of the re- 
maining questions concerns the behavior of slow waves, 
which were not present in the set up studied in this pa- 
per. In a high beta plasma (in which the sound speed 
is significantly larger than the Alfven speed) slow waves 
bear some similarities with Alfven waves, thus they may 
be expected to behave in a similar way. As their propaga- 
tion speed is somewhat slower than the Alfven speed, the 
surface where they would accumulate should be some- 
what beneath the Alfven surface. Because their group 
and phase velocity depend on their transverse structure, 
we could expect the region where they accumulate to be 
broader than in the case of Alfven waves. 

Based on the fact that the amplification of Alfven 
waves extracts energy from the stationary fiow, iWilliam"^ 
(fT975,) proposed that a broad turbulent zone would be 
created below the Alfven surface. Our analysis suggests 
instead that the amplification of Alfven waves energy flux 
creates an acoustic feedback propagating up and down 
rather than a turbulent zone. In addition to the inability 



of our one dimensional analysis to describe turbulence, 
the difference may be attributed to the fac t that we study 
a subsonic Alfven surface while iWilliami ()1975D studied 
a cold supersonic plasma, in which no upstream acoustic 
feedback is possible. 

Finally, we note that an Alfven surface could also ex- 
ist in an accretion fiow onto a magnetic object such as 
T-Tauri stars, magnetic white dwarfs or neutron stars. 
The dynamics may however be different in these con- 
texts especially if the Alfven surface lies in a supersonic 
flow. The possibility suggested by Williams (1975) that 
turbulence could play a role in these situations deserves 
further investigations. 
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APPENDIX 

WEAKLY COMPRESSIBLE ALFVEN WAVES 

In the linear regime, an Alfven wave perturbs the velocity and the magnetic field perpendicular to the stationary 
magnetic field as described in Equations (fT0| -(fT2 |) . On the other hand the perturbations of pressure, density, and 
longitudinal velocity and magnetic field are zero to first order. This is no longer true in the non linear regime where 
the Alfven wave becomes compressible. We estimate this with a second order development, where we take into account 
terms scaling like 5v^ and 5B^ but neglect all higher order terms. The second order perturbation of a variable X is 
denoted by 5^X, thus the background stationary flow is perturbed in the following manner : p p + 6'^p, v :^ v + S'^v 
etc... We then consider the conservation of magnetic flux, entropy, mass and energy at the interface between an 
unperturbed region and a region containing the Alfven wave. This interface is moving with respect to the observer 
frame at the propagation speed of the Alfven wave v + va + S'^{v + va)- The conservation of magnetic fiux imposes 
that the longitudinal magnetic field is unperturbed S^Bx — 0, and the conservation of entropy yields 6^P — p. 
The conservation of mass and energy across the interface can then be written : 

p&'^v-VAb'^p^^ (Al) 
- VA^''^ PVA + \vA {5Bl + SBD^Q iA2) 
Using these relations, one can then get the second order perturbations as a function of the linear wave amplitude 

S'P = c'6'p = - (A3) 
2/io (c^ - vj^) 



SB^ = 6Bl + 5Bl 



5^ 
2mo (c2 - v\) 
5'^p vaSB^ 



S v^VA = -- ^ 2T' (A5) 

THE LOW FREQUENCY REGIME 
Differential equation 

In a stationary state, the induction and the momentum equations are written (here, we use a resistivity and no 

viscosity however the same result would be found if one did the opposite assumption of neglecting the resistivity): 

= V X (v X B) +77AB, (Bl) 

VP (V X B) X B , ^ 

v.V v = B2 

P Pop 
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We linearize the problem by writing: 

v = w.Ux + (5v, (B3) 
B = S.Ux + (5B, (B4) 

where v and Bq are the velocity and magnetic field in the x direction in the unperturbed flow, while the perturbations 
and are assumed to be arbitrarily small. Assuming one dimensional perturbations (no dependence on y and 
z) and keeping only the first order perturbations, the y component of the induction and momentum equations can be 
written: 

which can be combined in one second order differential equation: 

dHBy ^ iv + v^)iv-v,) d6By ^ ^ ^ 

OX'' V ox ox 

In the vicinity of the Alfven surface, the profiles can be linearized: w ~ ug — 2io^x, va — vao + ^^vx and 2uj^x ^ vq. 
The differential equation then becomes: 

which is equivalent to Equation (|48p . The velocity perturbation can be obtained from the magnetic field perturbation 
as: 

c r> 

Svy = -VAo-^ (B9) 

Estimate of the maximum magnetic field 

In this subsection we estimate better the amplitude of the magnetic field accumulated at the Alfven surface without 
neglecting the linear coupling with a downward propagating Alfven wave. For this purpose we consider the two 
following conserved quantities: the transverse momentum, and the displacement of the field line. The downward 
propagating wave and the displacement accumulated at the Alfven surface should contain the same amount of these 
quantities as the incident Alfven wave: 

Sy-m = Sya + Syout (BIO) 

-PinWAin<52/in = -PoVAoSyo + PoutWAout^J/out (BH) 

Finally we get: 

s- Pout WAout + Pin ^^Ain f /T-)in\ 

dyo = ■ (B12) 

Pout^^Aout + PO^AO 

or equivalently: 

1 +WAout/l'Ain . miQ\ 

oyo = -;— -, oy-m (B13) 



1 + "WAout/wAO 

This result is compared successfully with the simulations in Figure 
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